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On  the  other  hand,  the  breadth  of  the  volurtK  distributions  for  the 
two  groups  did  not  increase  equally.  The  standard  deviation  of  the 
achiral  segment  volume  distribution  increased  by  1%  whereas  the 
chiral  group  distribution  widened  by  IS%.  It  can  be  seen  in 
Figures  2  and  3  that  the  largest  achiral  segments  found  are  about 
70  at  both  233’K  and  393'K.  On  the  other  hand,  the  maximum 
size  of  chiral  segments  increases  from  -70  A^  to  more  than  90  A 
The  smallest  segment  volumes  found  did  not  significantly  change 
for  either  type  of  segment.  ~ft  appears  that  raising  the  temperature 
induced  the  formation  of  a  fraction  of  chiral  segments  with  an 
exceptionally  large  volume  available  to  them. 

Flgvirt  3:  S«gm«nt  Volume  Distribution  et  393*K 


Table  I  presents  the  average  volume  of  the  chain  segments  which 
are  first  neighbors  of  achiral  and  chiral  se^ents  with  volumes 
that  pul  them  in  the  90tll  ,  9Slil,  and  991il  percentiles.  This 
information  gives  some  indication  of  how  these  larger  sized 
segments  are  distributed  in  space. 


233*K 

393*K 


Table  I 

Volume  PgKgntilc 


90 

achiral  chiral 
41.5  44.6 
43.7  45.6 


95 

achiral  chiral 

42.5  44.6 
43.9  45.8 


99 

achiral  chiral 
43.6  45.1 
44.5  46.1 


At  233’K  the  average  volume  of  the  neighbors  of  the  larger 
segments  was  somewhat  larger  than  the  average  for  segments  as  a 
whole  (39.0  A 7;  Figure  2).  This  indicates  that  the  segments  with 
the  largest  volumes  tended  to  be  associated  with  other  segments  of 
larger  than  average  volume. 

At  393'K  it  was  also  found  that  the  largest  segments  cluster  with 
segments  of  greater  than  average  volume  (42.3  A^;  Figure  3). 
However,  the  amount  by  which  the  average  volume  of  nearest 
neighbor  segments  was  enhanced  appears  to  be  less  than  that  at 
233'K.  This  implies  the  extent  to  which  clustering  occurs  is 
reduced  at  the  elevated  temperature. 

These  findings  agree  with  those  of  Deng  it  al.  [1 1]  Working  with 
an  idealized  two  dimensional  model  of  a  simple  atomic  glass,  they 
found  that  the  tendency  of  atoms  representing  an  excess  of  "free 
volume”  to  segregate  into  distinct  islands  was  greater  below  the 
glass  transition  than  in  the  melt. 

The  segmental  mobility  at  233'K  was  found  to  be  several  orders  of 
magnitude  lower  than  that  at  393'K.  [2|  The  observed  differences 
in  the  structure  at  the  two  temperatures  can  be  used  to  explain  this 
observation  qualitatively. 


The  calculated  Van  der  Waals  (O'K)  volume  of  the  achiral  CHl 
group  is  19.18  A^  while  that  of  the  chiral  CHR  group  is  38.21  A^. 
From  Figures  2  and  3  it  can  be  seen  that  virtually  all  of  the  achiral 
groups  have  more  than  their  Van  der  Waals  volume  associated  with 
them,  while  a  significant  portion  of  the  chiral  groups  have  less  than 
their  Van  der  Waals  volume  available  to  them.  Presumably,  this 
indicates  that  the  achiral  groups  are  relatively  free  to  move  while 
the  chiral  groups  are  more  constrained.  The  major  impediment  to 
large  scale  motion  of  the  polymer  backbone  therefore  seems  to 
the  chiral  groups. 

At  233'K,  22%  of  the  chiral  groups  have  less  than  their  Van  der 
Waals  volume.  These  groups  are  so  constrained  by  their  neighbors 
that  they  are  most  likely  unable  to  undergo  translations  or 
rotations.  The  average  length  of  chain  between  these  locked 
groups  is  -3  repeat  units.  It  has  been  postulated  on  the  basis  of 
mechanical  models  of  polymer  chains  and  other  considerations, 
that  the  smallest  length  of  unconstrained  chain  which  can  engage  in 
displacement  inducing  rearrangements  (i.e.  chain  mobility)  is  2-4 
repeat  units.  [9,  10]  If  this  is  so,  then  at  233'K  a  limited  number 
of  rearrangements,  which  can  yield  only  small  segment 
displacements,  are  available  to  the  unconstrained  lengths  of  the 
chain.  The  mobility  of  the  segments  will  therefore  be  low. 

At  393’K.  however,  only  10%  of  the  chiral  groups  have  less  than 
their  Van  der  Waals  volume.  The  average  length  of  chain  between 
locked  groups  is  consequently  -7  repeat  units.  The  number  of 
rearrangements  possible  and  size  of  the  displacements  they  produce 
are  increased  dramatically  and  hence  the  segmental  mobility  is 
likewise  enhanced. 
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Molecular  Dynamics  Simulation  of  Atactic 
Poly(propylcnc):  Structural  Difrerenccs  Between  the 
Liquid  and  Glass  State 

Mark  F.  Sylvester.  Sidney  Yip.  and  Ali.  S.  Argon 
Massachusetts  Institute  of  Technology 
Cambridge,  Massachusetts  02139 

IniCnuluctifln 

The  glass  transition  is  perhaps  the  single  most  profound  event 
occurring  in  polymers.  Despite  this,  it  is  still  poorly  characterized 
in  ail  but  a  phenomenological  manner.  In  particular, 
understanding  of  how  the  structure  of  a  polymer  glass  differs  from 
that  of  the  liquid  is  limited  to  the  somewhat  vague  concept  of  “free 
volume.”  Ill  It  IS  not  known  how  the  volume  associated  with 
atoms  might  be  distributed  within  the  polymer  or  how  changes  in 
the  distribution  are  related  to  the  property  changes  seen  at  the  glass 
transition.  In  addition,  what  other  structural  changes  occur  on 
passing  between  the  liquid  and  glass  states  are  not  known. 

We  have  previously  reported  on  a  molecular  dynamics  model  of 
amorphous  atactic  poly(propylene)  at  constant  temperature  and 
pressure/stress.  [2]  Initial  results  showed  that  a  variety  of  experi- 
mentaJ  thermodynamic,  structural,  and  transport  propenies  of  the 
polymer  in  both  the  liquid  and  glass  states  could  be  obtained  with 
reasonable  accuracy.  We  now  describe  the  use  of  the  model  in 
studying  the  detailed  microstructure  of  the  polymer  liquid  and 
glass. 

Metlunifllflgjt 

The  polymer  model  is  based  on  the  static  model  of  atactic 
poly(propylene)  developed  by  Theodorou  and  Suter.  [3]  It  consists 
of  a  single  polymer  chain  of  degree  of  polymerization  76  packed 
in  a  simulation  cell  with  periodic  boundary  conditions.  Hydrogen 
atoms  along  the  backbone  are  treated  explicitly,  while  a  unined 
atom  approximation  is  made  for  the  pendant  methyl  group. 
Realistic  potentials  are  used  for  bond  stretching  and  bending, 
intrinsic  barriers  to  rotation  about  backbone  bonds,  and  nonbondra 
dispersion  forces.  Parameters  for  the  potentials  are  based  on 
structural  and  spectroscopic  data  on  hydrocarbons.  Simulation 
under  constant  temperature  is  achieved  using  the  technique  of  Nosd 
[4],  while  constant  pressure  or  stress  simulation  is  carried  out  using 
the  methods  of  Andersen  [5)  or  Parrinello  and  Rahman  [6). 
respectively.  Further  details  of  the  simulation  may  be  found  in  the 
previous  paper.  [2! 

Initial  polymer  structures  above  (393'K)  and  below  (233'K)  the 
experimental  glass  transition  (253'K)  were  prepared  using  the 
energy  minimization  procedure  of  Theodorou  and  Suter  [3]  at  the 
experimental  polymer  density  for  the  two  temperatures  of  interest. 
The  atomic  coordinates  of  the  minimized  chains,  equivalent  to  a 
system  at  O'K.  were  used  as  input  to  the  molecular  dynamics 
simulation.  The  simulated  chains  were  brought  to  thermal  equilib¬ 
rium  at  either  233'K  or  393'K  by  first  carrying  out  a  10-15  ps 
microcanonical  simulation  with  occasional  velocity  rescaling  to 
bring  the  system  temperature  to  approximately  the  correct  value, 
followed  by  a  15-20  ps  canonical  simulation,  followed  by  a  70  ps 
isothermal-isobaric  simulation.  The  phase  space  trajectory  of  the 
chain  under  isothermal-isobaric  conditions  was  then  recorded  for 
20  ps  and  analyzed. 

The  micTOsmicture  of  the  polymer  was  probed  by  determining  the 
Voronoi  polyhedra  (7)  associated  with  each  atom.  The  polyhedra 
of  the  atoms  affiliated  with  each  chain  segment  (i  e.  CH2  or  CHR 
group)  were  joined  together.  The  distribution  of  the  volume 
contained  in  each  type  of  segment  polyhedra  (achiral  and  chiral) 
was  computed.  As  measure  of  the  tendency  for  large  polyhedra  to 
cluster,  the  volume  distribution  of  first  neighbors  of  those  poly¬ 


hedra  whose  volume  was  in  the  90111 . 9511l,  and  99111  percentiles 
were  evajuated.  In  addition,  the  distributions  of  bond  lengths, 
bond  angles  and  backbone  torsional  angles  were  calculated  from 
the  atomic  coordinates. 

Results  and  nUrM«if,n 

No  significant  differences  were  found  between  the  average  bond 
lengths  and  bonds  angles  in  the  polymer  structure  at  393'K  and 
233‘K.  However  the  distribution  of  bond  lengths  and  angles  was 
found  to  be  -30%  wider  at  393'K  than  at  233‘K.  This  is  consistent 
with  the  increased  motion  of  the  atoms  at  the  higher  temperature. 

Figure  1  show,  the  distribution  of  torsional  angles,  defined 
according  to  Flory's  convention  (81.  Substantial  differences  can  be 
seen  in  the  distribution  at  233'K  and  that  at  393'K.  An 
enhancement  of  the  fraction  of  torsional  angles  in  the  higher 
energy  g'  and  t*  states  occurs  at  the  expense  of  the  t  and  g*  states. 


Figure  1:  Distribution  of  Torsional  Angles 


The  distributions  of  Voronoi  polyhedra  for  achiral  and  chiral 
chain  segments  at  233'K  and  393'K  are  shown  in  Figures  2  and  3. 
respectively.  The  distributions  are  broad  and  asymmetric,  with 
pronounced  tails  at  large  volumes.  As  expected,  the  average 
volume  occupied  by  the  chiral  CHR  chain  segments  is  larger  than 
that  taken  up  by  the  achiral  CH2  groups. 

On  raising  the  temperature  from  233'K  to  393‘K,  there  were 
substantial  changes  in  the  dismbutions.  The  total  volume  of  the 
simulated  polymer  was  found  to  increase  by  8.6%.  The  average 
volume  of  the  achiral  segment  polyhedra  increased  by  9.3%.  while 
that  of  the  chiral  segment  polyhedra  was  raised  by  8.1%.  This 
slight  difference  is  probably  insignificant  and  indicates  that  both 
types  of  segments  contributed  equally  to  the  volume  expansion. 


0  10  20  30  40  30  00  70  OO  00  100 


Volum.  (A*  ) 


